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The general form of the electrostatic potential around an arbitrarily charged colloid at a flat 
interface between a dielectric and a screening phase (such as air and water, respectively) is analyzed 
in terms of a multipole expansion. The leading term is isotropic in the interfacial plane and varies 
with where d is the in-plane distance from the colloid. The effective interaction potential 
between two arbitrarily charged colloids is likewise isotropic and oc d~^, thus generalizing the dipole- 
dipole repulsion first found for point charges at water interfaces. Anisotropic, attractive interaction 
terms can arise only for higher powers d~" with n > 4. The relevance of these findings for recent 
experiments is discussed. 

PACS numbers: 82.70.Dd 



I. INTRODUCTION 

The self-assembly of stably trapped, sub-/im colloidal 
particles at water-air or water-oil interfaces has gained 
much interest in recent years. For the specific case of 
charge-stabilized colloids at interfaces, the repulsive part 
of their mutual interaction resembles a dipole-dipole in- 
teraction at large separations. This may be understood 
theoretically by approximating the colloid as a point 
charge located either right at the interface [l|, [2] (i-e. as- 
suming charges on the colloid- water interface) and/or 
above the interface Q (i.e. assuming charges on the 
colloid-air /oil interface). Additionally, the formation of 
metastable mesostructures with such colloids point to the 
possible existence of intercoUoidal attractions far beyond 
the range of van-der-Waals forces 0, however, care 
must be taken to avoid contaminations of the interface 
which lead to colloid mesostructures with similar appear- 
ance [t*]. Previous work % d, [T^, [HI, [ll aimed at re- 
lating this attractive minimum to capillary interactions 
due to interfacial deformations caused by a homogeneous 
surface charge on the colloids but with no conclusive an- 
swer. In recent work 5, 6], it was experimentally shown 
that the charge-carrying surface groups used for charge- 
stabilizing polystyrene colloids are actually distributed 
quite inhomogeneously and patchily over the colloid sur- 
face. Thus it was speculated in Refs. Q that through 
this inhomogeneous charge distribution like-charged col- 
loids could acquire effective dipole moments in the in- 
terface plane and attractive electrostatic interactions of 
dipole-dipole type could arise which might overcome the 
repulsion at shorter distances. 



* Supported by the German Science Foundation (DFG) via the 
Collaborative Research Center SFB— TR6 "Colloids in External 
Fields", project section D6. 



Motivated by the finding of inhomogeneous surface 
charge on colloids, we extend the asymptotic results 
for the electrostatic potential and interaction of point 
charges at water interfaces [3] to the general case of an 
arbitrary, localized colloidal charge distribution using a 
multipole expansion. The presence of the interface leads 
to restrictions in the multipole coefficients of the poten- 
tial around a single colloid and of the interaction energy 
between two colloids. In particular, we find that the lead- 
ing term in the effective interaction energy between two 
colloids at lateral distance d is isotropic in the interfa- 
cial plane, repulsive and oc regardless of the inhomo- 
geneities of the charge distribution in the colloids. Angu- 
lar dependencies enter the effective interaction potential 
only in higher orders. 



II. ELECTROSTATICS AT WATER 
INTERFACES 

A. A toy model: water as a perfect conductor 

For a quick insight on the effect of an interface on 
the multipole expansion of the electrostatic potential, 
we consider the water phase being a perfect conductor. 
The flat interface is located at z = and the colloid is 
modelled by a fixed charge distribution /9c(r) above the 
water phase. The boundary condition at z = simply 
implies that there is no tangential (or in-plane) electric 
field and the potential for z > can be obtained with the 
method of image charges. Therefore, the effective (real -I- 
image) charge distribution is spatially localized and can 
be enclosed in a ball of finite radius R (see Fig. [1] with 
— > 0). In standard spherical coordinates {s,9,ip) 
measured from the center of this ball, the potential in 
the upper phase for s > R can be written as a multipole 
expansion (in the remainder of the paper, the +(— ) index 
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FIG. 1: Geometrical configuration of the electrostatic prob- 
lem: the potential is calculated in the domain outside the 
sphere of radius R, which encloses the colloidal particle. The 
flat interface at z = separates an upper dielectric phase (di- 
electric constant e+) from a lower electrolytic phase (dielectric 
constant e-, Debye screening length k~^). 



will refer to evaluation in the upper (lower) phase): 



(1) 



in terms of normalized spherical harmonics Y^m- The 
boundary condition of vanishing in-plane electric field 
at the interface {0 = 7r/2) implies a^™ — for £ + m 
even. Thus, the monopole vanishes (ago = 0) as well as 
the in-plane dipole (ai±i = 0), and the leading decay 
is described generically by a dipole perpendicular to the 
interface (aio ^ 0). Consider a second, identical colloid 
located at an in-plane position d = (dx,dy). The total 
potential is now the linear superposition of the single- 
particle potentials by each colloid, and the electro- 
static energy of the two-particle configuration is 



U = U 







J d\pc{v) $*t(r + d) , 



(2) 



where U'^ is the energy in the limit d — > oo. Taylor 
expanding $° about r = one obtains to leading order 
in l/d 



(3) 



where Pz = 2 J d'^r z pc{r) = aio\/3/(47r) is the dipo- 
lar moment of pc and its image charge in the direction 
normal to the interface. This dipole-dipole interaction 
energy differs by a factor of one-half from the textbook 
result because 



B. Water as a conductor with linear screening 

The image charge construction in the case of perfectly 
conducting water provides an intuitive explanation of the 
origin of the normal dipole and the absence of an in-plane 
dipole. In the following we demonstrate that this find- 
ing still holds in the more realistic case of water being 
an electrolyte and the colloidal particle having an arbi- 
trary shape, possibly protruding into the region z < 0, 



with given charge distribution and dielectric properties. 
Assuming linear screening, the electrostatic potential sat- 
isfies (A — Kj.)^±{r) = with k+ = and — k being 
the inverse screening length in water. Using standard 
cylindrical coordinates (r, z,(p), we search for a solution 
outside a ball of radius R whose center is the coordinate 
origin and which encloses the colloid (see Fig. [T]) with 
the boundary conditions that the potential (i) vanishes 
at infinity, (ii) reduces to a given potential ^ii{9,ip) at 
the surface of the ball s = R, (iii) is continuous at the 
interface z = 0, and (iv) that the associated electric dis- 
placement perpendicular to the interface is continuous. 



I.e. 



9$^ 



dz 



dz 



{r > R). (4) 



The function ^ii{6^ip) is determined by the solution of 
the electrostatic problem inside the ball and contains the 
relevant information on the precise geometrical and elec- 
tric properties of the particle. By decomposing the prob- 
lem in the full domain into the solution of problems in 
simpler domains (the exterior of the sphere s = R and 
each of the halfspaces defined by z = (details can be 
found in App. E|, one can finally write the solution as 



the superposition $±(r) = ^'^^ {r, z,ip) -f <I>!^ {s,0,(p), 
where the contribution ^^''(r, z, ip) (using cylindrical co- 
ordinates) is given by 



+ 00 „oo 

$f(r,z,<^) = ^"""Z dqA,r,{q)J\„A{<ir)e 

m— — oo ^ 



(5) 

with K± — ± y^g^ -I- K j_ , and the contribution (us- 
ing spherical coordinates) reads 

oo +1 



=0 m=-i 



{sdsY 



(6) 



The coefficients are given by 



2tt 



±1 



$/? (61, 93) - $5''(r = i? sin 61, z = i? cos 61, ¥>) , (7) 



such that ^iP'' = at z = 0, $±(s ^ R) = ^R.iO,ip), 
and the boundary conditions (i)-(iii) are satisfied auto- 
matically. The coefficients A^{q) in the expression for 
<i>±'' (Eq. ([5])) must be chosen to enforce the boundary 
condition (|4]). This condition can be extended to the re- 
gion < r < i? by continuing the fields $±(r) with any 
virtual solution into the interior of the ball, s < R. The 
precise form of the continuation is irrelevant, since the 
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solution outside the ball depends only on the potential 
at the surface of the ball, ^ji{9,(p) (Faraday's cage ef- 
fect). Thus, by using orthonormality and closure of the 
set of Bessel functions { J|m| ('?'")}gG[o,oo) in the domain 
< r < oo, Eq. ([4]) can be solved for the coefficients 
Amiq)- 

Am{q) = ^ , , , ^, (8) 



with 

„± 



_e-'^^deY,,n{9 = TT/2,^)Ci^ and 

lirnil) Ir°° '^'^ T^fir) J\m\{qr), which are the Han- 
kel transforms of the radial dependence of the spherical 
part (see Eq. ^) continued into the region s < R 
by zero. Eq. ([5]) is not the explicit expression for the 
coefficients Am{q) because they appear implicitly also in 
the coefficients , see Eq. (O , but it does provide their 
dependence on q. In particular, for i — m odd (i.e., when 
^trn 7^ 0), the functions ■^f^{q) possess a Taylor expan- 
sion around q = Q with the lowest term being of order 
so that 

OO 

A,n{q) = ^ aj„iq\ with a-j^ = if j < |m|. (9) 

The existence of a Taylor expansion in q of the coeffi- 
cients Am{q) allows to extract the large-r behavior of 
the potential and the z-component of the electric field at 
the interface. Introducing the factors 

.CO _ 2^ r [ H+2±i A 
J,-™ lim / dppV|„,|(p)e-^P = /, , 

and inserting the expansion Q into the corresponding 
definitions of the fields, one obtains [l^l 



oo j-l 

$±(r,(^,z = 0) - ^ ^ ^ajmJ,m, (10) 
i=0 m=-j+l 

020 3 



dz 



aio -I- Cit 3 



where we have used that = if |to| = and C'^'^ = 
whenever Jj+i,„i — 0. Therefore, both the potential 
and the normal component of the electric field at the 
interface are asymptotically dominated by an angular- 
independent decay oc 1/r'^; anisotropic behavior arises 
only in subleading terms. By continuity, this conclusion 
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FIG. 2: Potential along the interface for an inhomogeneously 
charged cylinder half-immersed in water. Maxwell's equa- 
tions hold in the upper phase and the cylinder, the Poisson- 
Boltzmann equation holds in the lower phase. The parameters 
are given in the inset, the numerical calculations have been 
done using FEMLAB. 



also holds asymptotically for the fields at a fixed height 
h above or below the interface (r ^ \h\). 

This result is not exclusive of the single-particle con- 
figuration: if there are several particles at the interface, 
one can surround each of them by a ball of radius R and 
the solution <i>(r) of the electrostatic problem will be 
written as a superposition of single-particle potentials 
determined by the total potential at the surface of each 
ball (in general different from the potential $i? (0, (p) 
in the single-particle configuration). For each of these 
single-particle potentials the expansion ^ still holds, 
since it does not depend on the precise value of the 
potential at the balls. 



C. An illustrative 2d example 

We calculated the electrostatic potential for an 
inhomogeneously charged cylinder at an air-water 
interface (see the inset of Fig. [2] for some defi- 
nitions). Because of its two-dimensional nature, 
this problem is amenable to a numerical treatment. 
Here, the multipole expansion at the interface gives 
$±(z = 0) aoln|a;| -|- p^/x + q^x/x"^ + ■■■ and 

9$+/9zl^^Q ~ Pz/x^ -\- The numerical solution 

for (f>± shows that the in-plane dipole term oc 
is absent and the asymptotic expansion starts with 
the quadrupolar term (see Fig. [2]). The asymptotics 
for d^+/dz\^^Q (not shown) also contains the term 
oc which is interpreted as the effect of a counter- 
ion generated dipole component Pz perpendicular 
to the interface. These findings, most notably the ab- 
sence of Px^ match the previous ones in three dimensions. 
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D. The effective interaction energy 

The free energy functional of a multiparticle configu- 
ration in the linear screening approximation reads [l3| 



,(11) 



where the charge density pci^) is localized on the col- 
loidal particles. This includes the electrostatic energy 
as well as the entropy associated to the ion distribution. 
The extremum of F[^] provides the field equation in ther- 
mal equifibrium, V- [eV$] — en'^^ — Anpc- With the help 
of this equation, the free energy in equilibrium simplifies 
to 



^cq({x„}) 



dVpc(r)$(r), 



(12) 



which is known as the "potential of mean force" for the 
degrees of freedom (position of the center of a ball 
of radius R enclosing the a-th particle). One may de- 
compose Feq — F'^ + 6F, where F° is the equilibrium 
free energy in the limit |xq — X/jj — > oo (isolated par- 
ticles). The total potential can be similarly written 
as $ = J2a^a + where $°(r) denotes the poten- 
tial field generated by the a-th particle in isolation and 
i5$(r; {xa}) is the total perturbation induced by the pres- 
ence of other particles. Due to the linear nature of the 
problem, the perturbation (5$, 



(5$(r) = 



,|<-R 



G„^(r,r')$^(r'), 



(13) 



can be written in terms of a generalized susceptibility 
Gapiv, r') depending on the precise shape and charge dis- 
tribution of the particles. 

Since $°(r) near the interface exhibits asymptotically 
an isotropic decay cx 1/r^, (5<I>(r; {x^}) depends only on 
dap = |Xq, — x^l (and not on the orientation of Xq, — x^) 
in the asymptotic limit dap — *■ oo. Furthermore, (5$ is 
rescaled by a factor A~'^ if all distances dap are rescaled 
simultaneously by a factor A. From Eq. (fT^ the same 
property holds for 5F{{x.a})- In particular, for a two- 
particle configuration this yields an asymptotic potential 
of mean force of the form 



i^cq(d) - F° 



(14) 



and the constant B is positive for like particles. In anal- 
ogy with Eq. ([3]) , it is natural to interpret this expression 
as the interaction energy between two effective dipoles 
perpendicular to the interface. 

III. DISCUSSION AND CONCLUSION 

We have shown that the form of the multipole expan- 
sion of the potential around a charged colloid and of the 



effective interaction energy between two colloids trapped 
at a water interface is qualitatively different from the sit- 
uation in bulk. The dominating interaction terms can 
be qualitatively understood by assuming water to be 
a perfect conductor. The leading interaction term be- 
tween the colloids a distance d apart is of dipole-dipole 
type (oc d~^) and isotropic in the interfacial plane. In 
other words, even if the charges on the colloid surface 
are distributed arbitrarily the counterions arrange them- 
selves such that asymptotically the configuration corre- 
sponds to an effective dipole perpendicular to the inter- 
face. Orientation-dependent interactions and thus pos- 
sible attractions for like-charged colloids only arise in 
subleading order. 

This is in marked contrast to the analysis of the exper- 
iment reported in Refs. [1, @|. Motivated by the experi- 
mentally found inhomogenous surface charge, it was pic- 
torially suggested (see Fig. 1 in Ref. ^ ) that spontaneous 
fluctuations in the colloid's orientation would generate 
(via an instantaneously equilibrating counterion cloud) 
effective in-plane dipoles p.; with corresponding interac- 
tions oc [(i^(pi • P2) — 3(d ■ Pi)(d • p2)]/d^. After averag- 
ing over the orientation fluctuations, such an interaction 
would lead to an effective, isotropic attraction competing 
with the isotropic dipole-dipole repulsion. According to 
the model worked out by the authors, the total interac- 
tion potential would exhibit an attractive minimum due 
to the effect of the fluctuating in-plane dipoles at rather 
small distances [d ~ 2.2 colloid radii Rq, so small that 
already the use of a pure dipole-dipole interaction casts 
serious doubts on the reliability of the model). The anal- 
ysis in Ref. purported to support this picture is actu- 
ally incomplete and just states that no monopolar term 
arises, without entering into a systematic analysis of con- 
straints on higher order multipoles. In another note [Tsj 
the existence of the Taylor expansion of the coefflcients 
Am{q) around g = (see Eq. ©) was doubted on which 
the asymptotic analysis of the electrostatic potential and 
field is based. The present explicit proof of the analytic- 
ity of the coefficients Am{q) should disperse such doubts. 

The results of our work imply that asymptotically an 
in-plane dipolar interaction cannot arise if the counte- 
rions are equilibrated (see Eq. (|14p ). Consequently one 
cannot expect asymptotically relevant attractions from 
the orientational fluctuations of the colloids. However, 
for smaller d the asymptotic 1/d expansion is likely to 
break down. For small colloid radius, i?c ^ ^ this 
becomes relevant when d ^ in this case the screen- 
ing clouds of the colloids overlap and the interaction falls 
off exponentially with d before crossing over to the alge- 



braic decay 0, [Si- For large colloid radius, Rc S> k 
the precise shape and charge distribution of the colloids 
will determine the interaction whenever d ^ Rc- Cer- 
tainly, for both regimes a more elaborate numerical anal- 
ysis of the anisotropy in the electrostatic interactions is 
required to assess whether fluctuations in the orientation 
of the colloids may lead to attractions. However, even in 
that case their existence is doubtful looking at the gen- 
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eral results on the absence of like-charge attraction in 
confined geometries [l^. In any case, the results from 
the model studied in Refs. are not reliable since the 
model presupposes an interaction energy which does not 
satisfy the correct asymptotic decay given by Eq. (fT4|) . 



APPENDIX A: SOLUTION OF THE 
ELECTROSTATIC PROBLEM 

We consider the potential $±(r) in the domain shown 
in Fig. [1] given as the solution to the following problem : 



$+(r, (p,z ^0) 



0, r e {s > R, z > 0}, 

r e {s > R,z <0}, 

TT 



o<e<-, 



d2 



(f>_(r, (p,z ^ 0), 



z=0 

$(r)| < oo, 



dz 



< < TT, 

r>R, 
r> R, 



z=0 

r| — > oo. 



Here, ^ji{9,(p) is the potential at the surface of the ball 
s — R and is assumed to be given. In order to solve this 
problem, we split it in two auxiliary problems, one for 
each halfspace: 




z=0 



FIG. 3: Domain of definition of the problem UP. 



Problem UP in the domain s > R and z > 0, see Fig.|31 



0, 



r e {s > R,z > 0}, 



$+(r,v5,z = 0) = Fo{r,ip), r > R, 
|$+(r)| < oo, |r| ^ oo. 




z=0 



FIG. 4: Domain of definition of the problem LOW. 



Problem LOW in the domain s > R and z < 0, see 
Fig.m 



$_(.s = R,e,ip) 

$_(r,(^,z = 0) = Fo{r,Lp), 
$_(r)| < oo, |r 



r G {s > i?, 2 < 0}, 



2 

r > R, 
> oo. 



Here Fq (r, (p) (^potential at the interface 2; = 0) is an 
auxiliary function which will be eventually determined 
by the boundary condition ([4]). Each of these problems 
can in turn be decomposed in simpler problems, one 
with boundary conditions imposed only at the plane 
2 = and one with boundary conditions imposed only 
at the ball s = R: 



;Z=0 



FIG. 5: Domain of definition of the problem UP-cyl. 



Problem UP-cyl in the domain 2 > 0, see Fig. [SJ 



VHf = 0, 



r e {2 > 0}, 



$7'(r,(^,2 = 0) = Fo{r,ip)+To{r,ip), < r, 
|$7'(r)| < 00, |r| ^00. 



Problem UP-sph in the domain s > R, see Fig. (6] 



z=0 



FIG. 6: Domain of definition of the problem UP-sph. 




y2^sph 



0, re {s > i?}, 

= ^R{e,^)-<ff{R,e,^), o<e< 

p,sph 



TT 

2' 



|$;P"(r)| < 00, 



00. 



Here the function J-q (r, ip) verifies J-'q {r > R, tp) — and 
is an otherwise arbitrary smooth function which contin- 
ues the potential Fq (r, ip) into the region r < i? of the 
plane 2 = 0. As discussed in the main text, J-o{r,p) 
is just an intermediary auxiliary function whose precise 
choice is ultimately irrelevant for the determination of 
the total potential outside the ball s = R. With the 
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choice of boundary condition for $^ at s = i? it is clear 
that <i)^'^ = at z = and therefore 

$+(r) = $7' + $;P'', if re{s>i?,z>0}. 

Analogously, the problem LOW can be decomposed into 
a problem LOW-cyl and a problem LOW-sph and 

$_(r) = $1''' + $1P'\ if re{s>R,z< 0}. 

Each of these simpler problems is now amenable to an 



analytical solution, provided by Eqs. ©-([T]). The aux- 
iliary functions Fq and J-q are absorbed in the unknown 
coefficients A,m{q) in Eq. ([5]), which are then determined 
by Eq. the only boundary condition of the original 
problem not taken into account by the stepwise process 
of decomposing the problem into simpler ones. 
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